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Perturbation Theory on the Superconductivity of Heavy Fermion 
Superconductors CeIr x .Coi_. r In5 
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We reformulate the Eliashberg's equation for the superconducting transition within the 
quasi-particle description, which takes account of the heavy electron mass. We discuss the 
superconductivity of Celr^Coi^^Ins by using such a renormalized formula. On the basis 
of an effective two-dimensional Hubbard model which represents the most heavy quasi-two 
dimensional Fermi surface with /-character of Celr^Coi-^Ins, both normal and anomalous 
self-energies are calculated up to third order with respect to the renormalized on-site re- 
pulsion U between quasi-particles. The superconducting transition temperature is obtained 
by solving the Eliashberg's equation. Reasonable transition temperatures are obtained for 
moderately large U . It is found that the momentum and frequency dependence of spin fluc- 
tuations given by RPA-like terms gives rise to the d-wave pairing state, while the vertex 
correction terms are important for obtaining reasonable transition temperatures and for dif- 
ferentiating n-dependence of T c . These obtained results seem to explain the x-dependence 
of T c in Celr^Coi-jJns and confirm the usefulness of the perturbation theory with respect 
to U. 

KEYWORDS: heavy fermion systems, electron correlation, Ce(Ir, Co, Rh)In5, superconducting 
transition temperature, perturbation approach, vertex correction 

1. Introduction 

In heavy fermion system, the effective mass is enhanced by 10 2 ~ 10 3 times as large as that 
of free electron. The coefficient of T 2 -law of resistivity is also enhanced by 10 4 ~ 10 6 times 
as large as that of conventional metal. Such a large electron mass and coefficient of T 2 -law 
of resistivity stem from the electron correlation. In heavy fermion superconductors, itinerant 
electrons compose such a Fermi liquid and then undergo the superconducting transition. 
From this point of view, we can expect that superconductivities in heavy fermion systems are 
derived also from the electron correlation through the momentum and frequency dependence 
of the effective interaction between electrons. Therefore, in contrast to the cuprates, it is not 
clear whether taking only the effect of the specific diagrams which are favorable to the spin 
fluctuations is reliable for revealing the mechanism of superconductivity in the heavy fermion 
sup er conductor s . 

In this paper, we present an essentially important point in theoretical study of supercon- 
ductivity in heavy fermion systems. In contrast to the cuprates, the momentum dependence 
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of the mass enhancement factor in heavy fermion systems is weak because of the localized 
nature of /-electrons. In this case we can reformulate the Eliashberg's equation also in the 
renormalized form, where we have already included the large mass enhancement effect. By us- 
ing the renormalized Fermi energy as unit of energy the critical temperature T c determined by 
the perturbation is the same as that obtained after taking account of the mass enhancement 
effect. Thus, we can formulate the perturbation theory which is consistent with the large mass 
enhancement. The formulation is one of purposes of the present paper. Therefore, we calcu- 
late the superconducting transition temperature T c of heavy fermion superconductors by the 
perturbation theory based on Fermi liquid theory. In our formalism, the effective interaction 
corresponding to the superexchange interaction 1 is included in the vertex terms in principle, 
because we start from renormalized Fermi liquid and treat only the Coulomb repulsion by the 
perturbation theory. 

The perturbation approach is sensitive to the dispersion of the bare energy band by its 
nature, it implies that the lattice structures and the band filling play the essential roles in 
the calculation of T c . Therefore it is important to evaluate T c on the basis of the detailed 
electronic structure in each system. In this paper, we calculate T c of Celr^Coi-^Ins. 

The organization of this paper is as follows. In §2, a theoretical treatment of heavy fermion 
superconductivity is presented. In §3, T c of Celr^Coi-^Ins is calculated, after giving the 
model Hamiltonian, perturbation expansion terms and Eliashberg's equation. Moreover the 
selfenergy, density of states and other quantities are shown. Finally in §4, summary, discussions 
and conclusion are presented. 

2. A theory of heavy fermion superconductivity 

In this section, we discuss a theoretical treatment of heavy fermion superconductivity. As 
indicated in various experiments, such as a jump of the specific heat at the transition tem- 
perature (T c ), the superconductivity can be interpreted as a transition of quasi-particles with 
heavy electron mass. We here formulate the way to calculate the superconducting transition 
temperature T c on the heavy fermion quasi-particle state in the periodic Anderson model 
(PAM). 

One of the characteristics in heavy fermion systems is the weak momentum dependence, 
corresponding to degree of the local /-electron spins at high temperatures or energies. The 
perturbation expansion in the hybridization matrix element in PAM leads to the RKKY 
interaction between these local spins. This expansion is valid at high temperatures, and may be 
the origin of the magnetic transition in /-electron systems. However, if no magnetic transition 
occurs, the system goes to a singlet ground state as the whole. This just is the Fermi liquid state 
with heavy electron mass. In this case, we had better treat it by the perturbation theory with 
respect to the on-site Coulomb repulsion U between /-electrons, because of the analyticity 
about U. The analytic property has been exactly proved in impurity Anderson model. 2-5 
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Although there is no exact such proof in PAM, the principle of the adiabatic continuation 
confirms the correctness of the perturbation theory in U. 

From such a point of view, one of the authors (K.Y.) has described the Fermi liquid theory 
for the heavy fermion state in PAM. 6 The momentum independent part of the 4-point vertex 
functions r)°S, between /-electrons with the spin a and a' plays an important role. This s- 
wave scattering part, corresponding to the above-mentioned local character at high energies, 
leads to the featureless large mass enhancement factor 7 = a -1 . In fact, when the momentum 
dependence of the mass enhancement factor can be ignored, the T-linear coefficient of the 
specific heat in PAM, is given by 7 ~ (2/S)ir 2 k^p f (0) 2 F l ff = (2 / '3)vr 2 k\p j (0)7, where k B is 
the Boltzmann's constant, and p* (0) the /-electron density of states at the Fermi level. This 
indicates that the large mass enhancement factor 7 = a -1 is represented as pi* (0)r 1 T ° J c with 
use of the s-wave scattering part of the 4-point vertex. Thus, we can see that the vertex F l °£ 
is enhanced by 7 = a -1 . This corresponds to the fact that the interaction a 2 F l °£ between 
the quasi-particles has the order of magnitude of the effective band-width of quasi-particles 
To ~ 1/7/9-' (0); To is the characteristic energy scale of the heavy fermion state, and below To, 
the low energy excitation can be described by the Fermi liquid theory. On the other hand, 
the imaginary part of the /-electron self-energy, which is proportional to the T 2 -term of the 
electrical resistivity, is given by A fc = (47r 2 /3)p / (0) 2 7rr 1 T ° J c T 2 , if the momentum depend ence of 
the vertices can be ignored. Thus, the coefficient A of the T 2 -term of the electrical resistivity 
is proportional to 7 2 through the vertex rj° c - This relation A oc 7 2 is just the Kadowaki- 
Wood's relation, which holds in the case where the momentum dependence of the vertices is 
sufficiently weak to be ignored. Actually, in the heavy fermion systems, the Kadowaki- Wood's 
relation is confirmed. This fact means that the 4-point vertex function, i.e., the interaction 
between quasi-particles possesses the above-mentioned local nature. 

The superconductivity in the heavy fermion systems appears under this situation below 
To. The s-wave large repulsive part rl^ prevents the appearance of the s-wave singlet super- 
conductivity in strongly correlated systems. In this case, an anisotropic (such as p-, d-wave 
and so on) superconductivity will be realized due to the remaining momentum dependence of 
the 4-point vertices. We here discuss such a transition to the anisotropic superconductivity on 
the basis of the heavy fermion quasi-particle state renormalized by the large s-wave scattering 
part 

First of all, we divide the 4-point full vertex function T crcr >(pi,p2;p3,P4) into the large 
s-wave scattering part F l ™, and the non-s wave part Ar CTCT /(pi,p2;P3,P4); 

^aa'(Pl,P2-,P3,P4) = F 1 ™, + AY afT i ,P2'-,P3: Pi)i (1) 

where pi andj<2 are the incident momenta, and P3, p^ the outgoing. The momentum dependent 
part has rather remarkable momentum dependence as the heavy fermion quasi-particles grow 
below Tq. The effective Cooper pairing interaction in the anisotropic superconductivity results 
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from the momentum dependence of AT aa i(pi, P2]P3,P4)- Our purpose is to formulate how to 
calculate AT <J<J i(pi,p2;p3,P4) for the heavy fermion quasi-particles, which are renormalized 
by the large local part 1^°°- Corresponding to these two terms of the full vertex function, we 
can also separate the self-energy into the local part and the non-local part 

E(k,u) = S 1oc (w) + AS(A;,lj). (2) 

In this case, the /-electron Green's function below To is given by 

1 



G(k,u) = 



oj-E k - t(k,u) - 



+ GincH, (3) 



u>-e k 

where a is a wave function renormalization factor, and = a(£fc+Re£ loc (0)), = \faVk and 
E(fe, to) = aAT,(k, u) are, respectively, a renormalized /-electron band, hybridization term and 
self-energy. G- mc {uj) is a featureless incoherent part. The imaginary part of S loc proportional to 
lo 2 can be included in the remaining renormalized self-energy T,(k,u>), as indicated by Hewson 
for the impurity Anderson model. 7 ' 8 From the Ward-Takahashi identity, we can obtain the 
large mass enhancement factor, 

7 = a" 1 = 1 - = l-lj ^Gl c (u). (4) 

The above equation of the Green's function is the counterpart on the quasi-particle description. 
In addition, we can set the large local vertex T l °£ as a constant value at zero frequencies, as far 
as the frequency of external lines is less than T . In this case, the effective on-site interaction 
a 2 r 1 °j C works between the quasi-particles. In the impurity Anderson model, a 2 T^ = an A = 
4Tk, where A and Tk are, respectively, the width of the virtual bound state and the Kondo 
temperature. Therefore, we can assume that a 2 T l °^ has the order of To. This is also consistent 
with the above-mentioned discussion of the coefficient of the specific heat. 

We now show that in the region u < To, the renormalized self-energy Y±(k,u>) and vertex 
function T a a'(Pi,P2',P3iP4) = o? Ar ' a a'{pi-,P2',P3,Pi) can De discussed with the perturbation 
expansion with respect to the effective on-site interaction F^ = a 2 T l °£. This expansion can 
be also considered as the one with respect to the momentum dependence which is governed by 
the renormalized Fermi surface. First, let us consider the renormalized self-energy Y±(k,u>) = 
aAY,(k,oj). Because we have divided the Green's function into the two parts aG(k,uj) + 
Gi QC (u>), we can correspondingly divide any diagrams in the expansion with respect to the 
bare interaction U into diagrams with and without the momentum dependence. The latter is 
contained in the local self-energy Si oc (lj), on the other hand, the former can be rewritten as 
the expansion with respect to the effective interaction r 1 ^ by concentrating diagrams with 
the same momentum dependence. As shown in Fig. 1, for instance, the C/ 2 -type diagram has 
the same momentum dependence as the T^-type, each vertex of which includes all the s- 
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wave scattering process. Because three lines of the renormalized Green's function aG(k, to) 
yields a 3 , and the renormalized self-energy T,(k,u) includes a, the ri^-type is estimated by 
calculating a^T^^GGG. This just corresponds to the diagram of the order two with respect 
to the effective interaction F^P? = a?T]P?. In the same way, the full vertex T acr i(pi,P2',P3,P4} 
can be also re-written as the expansion with respect to the quasi-particle Green's function 
G{k, lo) and the quasi-particle s-wave scattering interaction r)°S, as shown in Fig. 2. Thus, 
we can discuss the phenomena below To within the renormalized effective PAM, so far as the 
frequency of external lines is less than To. The validity of such a perturbation expansion with 
respect to the renormalized interaction is shown by the fact that the fixed point Hamiltonian in 
the renormalization group can be written as the renormalized PAM. In this case, the locality of 
the relevant interaction is very important. In the heavy fermion systems, this condition almost 
holds as shown in the Kadowaki- Wood's relation. We formulate below the superconducting 
transition on the quasi-particle description discussed above. 

The superconducting transition is always marked by divergence of the full vertex 
^aa'(p, ~P',p' i ~p') = ^aa'(p,p')- We can rewrite it as divergence of the quasi-particle full 
vertex F aa i(p,p'). As usual, we estimate the Eliashberg equation, 

J2r {2) (P,p')\G(p',iu;' n )\ 2 A(p',iuj' n ) = A(p,iu n ), (5) 
P > 

where A(p, ito n ) is an anomalous self-energy with a fermion Matsubara frequency io n = 
(2n + l)vrT, and the Cooper pairing effective interaction T^ 2 \p, p') is the particle-particle 
irreducible vertex. The Eq.( 5 ) includes the integral of \G(p', iuj' n )\ 2 . The most important 
part of this integral comes from the part mediated by quasi-particles a 2 \G(p', iuj' n )\ 2 . The 
Eliashberg equation can be re-written as 

Y j a 2 T^{p,p')\G{p',iu' n )\ 2 A{p\iJ n ) = A(p,iu n ). (6) 

This a 2 rW{p ,p') is the particle-particle irreducible vertex between quasi-particles with the 
external frequencies less than To. As discussed above, for such a vertex, we can also apply 
the perturbation expansion with respect to the renormalized on-site repulsion a 2 r!j.° c between 
quasi-particles. 

In order to treat the heavy fermion superconductivity, we have approximately introduced 
the quasi-particle state renormalized by the s-wave scattering part, which itself does not 
yield the pairing interaction. Starting from the quasi-particle state, we calculate the momen- 
tum dependent interaction between quasi-particles by the perturbation theory and derive the 
anisotropic superconductivity. It should be noted that even if the effective interaction depends 
on momentum as seen in the cuprates, we can treat such a separation of the vertices in the 
momentum space, as far as, at the final stage, everything is included in a consistent way. 
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3.1 Introduction to Ce(Rh, Ir, Co)In 5 

Superconductivities are discovered in a series of CeMIns tetragonal compounds for 
M=Co, Rh or Ir. Celrlns 9 and CeCoIns 10 are superconductors at ambient pressure. The tran- 
sition temperature T c (and the electronic specific heat coefficients 7) are 0.4K (680mJ/K 2 mol), 
and 2.3K (300-1000mJ/K 2 mol), respectively. In contrast, CeRhlns 11 orders antiferromagneti- 
cally below Tn = 3.8K, whereas the superconductivity is observed under pressure, p > 1.6GPa. 
The observations of power law temperature dependence of the specific heat, thermal conduc- 
tivity, and NMR relaxation rate have identified Ce(Rh, Ir, Co)In5 as unconventional supercon- 
ductors with line nodes. The thermal conductivity measurement in a magnetic field rotating 
within the 2D planes also reveals that the superconducting gap symmetry of CeCoIns most 
likely belongs to d x 2_ y 2. 12 The temperature-composition phase diagram for Ce(Rh, Ir, Co)In5 
has been obtained mainly from heat capacity measurements on single crystals. 13 The phase 
diagram shows an interesting asymmetry with respect to the coupling between magnetism and 
superconductivity: in the case of Rh-Ir, T c increases as the magnetic boundary is approached, 
whereas in Rh-Co case, T c decreases. According to the phase diagram, the superconducting 
phase of Celr^Coi-^Ins is apart from the antiferromagnetic ordered phase. When x increases 
from to 1, T c decreases from 2.3 K to 0.4 K. The temperature dependence of the resistivity 
p of Celrlns follows the dependence p = po + AT 2 . There exists a large anisotropy in the 
resistivity of Celrlns. This anisotropy is closely related to the quasi- two dimensional Fermi 
surface. The band calculations were performed for Celrlns 14 and CeCoIns , 15 The Fermi sur- 
faces calculated for Celrlns and CeCoIns, respectively are almost same. The topology of the 
Fermi surface is explained by the 4/-itinerant band model. The calculated Fermi surfaces 
consist of nearly cylindrical Fermi surfaces and small ellipsoidal ones. 

3.2 Model Hamiltonian for Celr^Coi-^Ins 

We calculate here T c of Celr^Coi-^Ins. According to the discussion in §2, we start from the 
quasi-particle state and then calculate the momentum dependent interaction between quasi- 
particles by the perturbation theory and derive the superconductivity of Celr^Coi-^Ins. Tak- 
ing the importance of the most heavy quasi-two dimensional Fermi surface with /-character 
into account in realizing the superconductivity in Celr^Coi-^Ins, we represent the band as 
an effective two-dimensional Hubbard model on the tetragonal lattice. In this section, we call 
quasi-particle "electron" and call the renormalized on-site repulsion between quasi-particles 
fioc " Coulomb repulsion U " . We consider only the nearest neighbor (effective) hopping inte- 
grals t. We rescale length, energy, temperature, time by a, t, | respectively (where a,kB,ft 
are the lattice constant of the tetragonal basal plane, Boltzmann constant, Planck constant 
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divided by 2ir respectively), we write our model Hamiltonian as follows: 

H = H + H 1 , (7) 

H = ^(e(k)-^)4 CT a k(J , (8) 

k,cr 

e(k) = 2(cos k x + cos k y ), (9) 

Hl= 2N^2 ^2 5 ki+k2,k3+k4 a ic lCT a k 2 a' a k3^ a k4 CT > ( 10 ) 
cr^cr' kik2k3k4 

where aj^Oko-) is the creation(annihilation) operator for the electron with momentum k 
and spin index a; t and \i are the hopping integral and the chemical potential, respectively. 
The sum over k indicates taking summation over a primitive cell of the inverse lattice. The 
Fermi surfaces calculated for Celrlns and CeCoIns, respectively are almost same. Therefore 
we consider that the electron number n per one spin site and the Coulomb repulsion U are 
determined corresponding to the composition x. 

3. 2. 1 model parameters 

Our model parameters are the Coulomb repulsion U and the electron number n per one 
spin site. The bandwidth Ae of our model is 8. According to the band calculation and the 
de Haas- van Alphen experiment, 14,15 the parameter region of our model which reproduces 
the Fermi sheet is given by 0.345 < n < 0.385. We consider that n = 0.345 and n = 0.385 
are the best fitting values of parameters which well reproduce the Fermi sheet of Celrlns 
and CeCoIns, respectively. The superconducting transition temperatures T c of Celrlns and 
CeCoIns determined by experiments are about 0.4 K and 2.3 K, respectively. To evaluate the 
rescaled T c (= k^Tc/t) , we have to estimate the renormalized bandwidth Ae(= 8t) because 
we have started from the renormalized Fermi liquid. We estimate Ae by using following equa- 
tion; Ae ~ ra^nd Aebandi where m-bandi m <mvA, Aeband are the effective mass calculated by the 
band calculation, 14 ' 15 the cyclotron effective mass obtained the de Haas- van Alphen by ex- 
periment, 14 ' 15 the bandwidth calculated by the band calculation, 14 ' 15 respectively. Then the 
rescaled superconducting transition temperatures of Celrlns and CeCoIns are, respectively, 
T c (CeL-in 5 ) — 3.2 x 10~ 3 and T c( - CeCoInr) ) ~ 3.3 x 10~ 2 when we set 0.1 Ry as the value of 
Ae band n5 (- Ae band° In5 ) according to the band calculation. 14 

3.3 Green's functions 

3. 3. 1 Bare Green's functions 

The bare Green's function is the following, 

G<°>(fc) = : ]- , (11) 
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where e n = (2n + l)nT(n : integer) is the fermion-Matsubara frequency and short notation 
k = (k, e n ) is adopted. We consider the Hartree term is included in the chemical potential. 



3.3.2 Dressed normal Green's functions 

Next we consider the dressed normal Green's function. When we consider the situation near 
the superconducting transition temperature, the Dyson-Gorkov's equation can be linearized. 
Therefore the dressed normal Green's function is obtained from the bare Green's function 
with only the normal self-energy correction. We expand the normal self-energy up to third 
order with respect to J7;the diagrams are shown in Fig. 3. 

Then we obtain the normal self-energy as follows, 

= ^ ZKxot* " + U ^ k ~ k "> + U3 ^ k + WPotf), ( 12 ) 
k' 

where xo( - ■ ■ ) an d (po(- • • ) are given respectively as 

Xo(q,uj m ) = -— ^G (k, e n )G (q+k,cj m + e n ), (13) 

k,n 

0o (q, w m ) = -— ^2 G (k, e n )G (q - k, u> m - e n ). (14) 

k,n 

Here u: m = 2mnT{m : integer) is the boson-Matsubara frequency. The quantity xo^^m) has 
the physical meaning of the bare susceptibility and expresses spin fluctuations in the system. 
More over, the bare susceptibility plays an important role in the calculation of T c , that is, it 
determines the magnitude and the spatial and temporal variation of the effective interaction 
between electrons, through the higher order terms in U. (See the first term of right hand side 
of the equation 20.) 

Note that the Hartree term has been already included in the chemical potential and the 
constant terms which have not been included in the Hartree term are included in the chemical 
potential shift when we fix the particle number. 

Then the dressed normal Green's function is 

G(k) = ie n -(e(k) -//-<*,/ + £„(*))' (15) 
where d/j, is determined so that the following equation is satisfied. 

k k 

We expand the above equation up to the third order of the interaction with regard to S/j, — 
S n (fc), we obtained S/i as 

, M= _5E»y)s»w (17) 



Xo(0,0) 
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3.3.3 Anomalous self- energy and effective interaction 

When we consider the situation near the superconducting transition temperature, the 
anomalous self-energy T, a (k) is represented by the anomalous Green's function F(k) and 
the (normal)effective interaction. We expand the effective interaction up to third order with 
respect to U as shown in Fig. 4. 

Then we obtain the anomalous self-energy as follows; 

= -^^ U + U2 ^ k + k ') + 2U ^ k + k ')^ k ') 
k' 

~ U3 ^2 E G o(k')[xo(k + k') - Mk + k')]G {k + k' - k")F{k") 

k'k" 

~ U3 ^2 E G (k')[ X o(-k + k') - M-k + k')]G (-k + k! — k")F{k"). (18) 
k'k" 

3.4 Eliashberg's equation 

From the linearized Dyson-Gorkov equation, we obtain the anomalous Green's function 
as follows; 

F(k) = \G(k)\ 2 ^ a (k). (19) 
Then the Eliashberg's equation is given by 

E «( fc ) = -^E^ + ^o^ + ^ + ^o^ + ^IG 1 ^)! 2 ^^) 
k' 

- U 3 ^ J2 G o(k')[xo(k + k') - Mk + k')]G (k + k> - k")\G{k")\ 2 Y> a (k") 
k'k" 

~ U '^2 E G (k')[ X o(-k + k') - M-k + k')]G (-k + k> - k")\G{k")\ 2 V a {k"). 
k'k" 

(20) 

This equation is corresponding to the Eq.( 6 ). We consider that the system is superconducting 
state when the eigen- value of this equation is 1. 

3.5 Calculation Results 

3.5.1 Details of the numerical calculation 

To solve the Eliahberg's equation by using the power method algorithm, we have to cal- 
culate the summation over the momentum and the frequency space. Since all summations are 
in the convolution forms, we can carry out them by using the algorithm of the Fast Fourier 
Transformation. For the frequency, irrespective of the temperature, we have 1024 Matsubara 
frequencies. Therefore we calculate throughout in the temperature region T > Tu m , where 
Tii m is the lower limit temperature for reliable numerical calculation, which is estimated about 
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2.0 x KT 3 (> Ae/(2vr x 1024) ~ 1.2 x 10~ 3 ); we divide a primitive cell into 128x128 meshes. 
We have carried out analytically continuing procedure by using Pade approximation. 



3.5.2 Dependence of T c on U,n and vertex correction terms 

To solve the Eliashberg's equation, we set the initial gap function {d x i_ y i -symmetry) as 
follows. 

T> a (k) oc cos k x — cos k y . (21) 

The calculated gap functions show the node at k x = k y and k x = —k y and changes the 
sign across the node in all approximations and for all parameters. The symmetry of Cooper 
pair is d x 2_ y 2. The dependence of T c on U, n are shown in Fig. 5. To examine how the ver- 
tex corrections influence T c , we also calculate T c by including only the RPA-like diagrams of 
anomalous self-energies up to third order, in other words, without the vertex corrections. We 
compare obtained T c (RPA-like only) with T c calculated by including full diagrams of anoma- 
lous self-energies up to third order, in Fig. 5. From this figure, we can point out the following 
facts. For large U higher T c are obtained commonly for all cases. T c calculated by including 
only the RPA-like diagrams of anomalous self-energies up to third order is higher than T c 
calculated by including full diagrams of them. These results show that the main origin of the 
d-wave superconductivity is the momentum and frequency dependence of spin fluctuations 
given by the RPA-like terms. The vertex corrections reduce T c by one order of magnitude. 
So the vertex corrections is important for obtaining reasonable T c and for differentiating n- 
dependence of T c . When we fix the Coulomb repulsion U and increase n(0.345 ~ 0.385), the 
system get close to the half-filling state (n = 0.5). In this case, the momentum dependence of 
Xo(Q)O) is slightly enhanced (see Fig. 6 shown later) and at the same time Fermi level gets 
close to the van Hove singularity (see Fig. 7 shown later), then higher T c are obtained. If we 
assume £/ceirin 5 — ^CcCoins, T c of Celrlns with n = 0.345 is lower than that of CeCoIns with 
n = 0.385 and this tendency is in good agreement with experimental results. 

3. 5. 3 Behavior of xo (q, 0) 

The calculated results of the static bare susceptibility are shown in Fig. 6 for various values 
of n. From this figure, we point out the following facts. In the case of half-filling state, xo(q> 0) 
has commensurate peak at (tt,tt). When we increase n(0.3 ~ 0.5), the system get close to the 
half-filling state and the peak and the momentum dependence of Xo(q> 0) are enhanced. In the 
case of n = 0.345 ~ 0.385, Xo(q>0) has the incommensurate peak around (vr,7r). This peak 
is not prominent but Xo(q>0) has sufficiently strong momentum dependence. These results 
mentioned above indicates the d x 2_ y 2 symmetry of the gap function. 
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3.5.4 Density of states 

The density of states(DOS) is given by 



k 



(22) 



where 



G K (k,u,) = G(k,e„)| ; 



We show the n-dependence of DOS in Fig. 7. From inset in this figure, we can see that the 
position of the van Hove singularity shifts upward from the Fermi level with decreasing n. 
This departure of the van Hove singularity from the Fermi level reduces the superconducting 
transition temperature T c . 

3.5.5 Self- energy 

The self-energy is given by 



The real part and the imaginary part of the self-energy at Fermi momentum are shown 
in Fig. 8 and Fig. 9 respectively. The w-dependence of both parts near to = are respectively 



and the damping rate of the quasi-particle become larger as U increases. These results are the 
typical Fermi liquid ones. 

4. Summary, Discussion and Conclusion 

In this paper, we have reformulated the superconducting transition on the quasi-particle 
description and discussed of the superconductivity of Celr^Coi-^Ins on the basis of such a 
renormalized formula. By the low order perturbation expansion it is difficult to obtain ac- 
tual mass enhancement factor. Here we have started from heavy fermions possessing a large 
effective mass. To obtain the anisotropic superconductivity we have taken account of the mo- 
mentum dependence of interaction between quasi-particles. By this procedure we can make an 
argument consistent with both heavy electron mass and anisotropic superconductivity. Taking 
the importance of the most heavy quasi-two dimensional Fermi surface into account in realiz- 
ing the superconductivity in CeIr. E Coi_. r In5, we have represented these bands with an effective 
single band two-dimensional Hubbard model and calculated T c on the basis of the third order 
perturbation theory. We have pointed out that the main origin of their superconductivities can 
be ascribed to the momentum and frequency dependence of the spin fluctuations, although 
the vertex correction terms are important for reducing T c then obtaining reasonable transition 
temperatures and differentiating n-dependence of T c . 




(k, W ) 



— S n (k, e n )|j 



given by Re£^(k^,u;) oc — oj and ImS^(kj,cj) oc —uj 2 . This behavior is the same as that for 
the usual Fermi liquid. As U increases, the slope of ReX^(ky,tj) at oj = becomes steeper 
and the coefficient of the u; 2 -term in ImS^(kj, u>) become larger. This indicates that the mass 
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Almost the same results have been obtained in the calculation of T c for high-T c cuprates 
performed by Hotta 16 and for organic superconductors performed by Jujo et al. , 17 based on 
the third order perturbation theory. For high-T c cuprates, Hotta calculated T c on the basis of 
d-p model and succeeded in reproducing experimentally obtained dependence of T c on carrier 
number for overdoped cuprate within the third order perturbation theory. In spin-triplet 
superconductor Sr2Ru04, Nomura and Yamada have recognized that the momentum and 
frequency dependence of the effective interaction between electrons which is not included in the 
interaction mediated by the spin fluctuation is essential for realizing the spin triplet pairing, by 
treating the electron correlation within the third order perturbation theory, on the basis of the 
two-dimensional single-band Hubbard model. 18 According to their theory, third order vertex 
correction terms, which are not direct contribution from spin fluctuations, are important for 
the superconductivity in Sr2Ru04. They also searched finite values of T c within the second 
order perturbation theory (the second order term is the only term which is included in the 
interaction mediated by the spin fluctuation within the third order perturbation theory in 
the case of spin-triplet superconductivity of the single-band Hubbard model) and the random 
phase approximation, but could not find any finite value of T c within the precision of their 
numerical calculations and for appropriate values of parameters. Recently Kondo has examined 
the superconductivity of the two-dimensional Hubbard model in the small U limit within the 
second perturbation theory. One of his conclusion is also that the spin-triplet superconducting 
state is not stable in the small U limit. 19 

These facts described above suggest that the calculations of T c which include only the 
spin fluctuation terms are questionable and should be carefully performed with other terms. 
In our calculation, the effects of other terms are taken into account by third order vertex 
correction terms which are not included in the interaction mediated by the spin fluctuation. 
We assume that the effect of the electron correlation is almost taken into account by the third 
order perturbation theory. To confirm this assumption, we have to do with the forth order 
perturbation theory. The forth order perturbation theory is a remaining future problem. We 
point out that in calculating the superconducting transition temperature T c of typical heavy 
fermion superconductors, it is reasonable to treat the electron correlation by perturbation 
theory based on the Fermi liquid theory instead of taking only the effect of the specific 
diagrams which are favorable to the spin fluctuation. 

The electronic structure of pressured CeRhIns(P > I.6GPa) is considered to be same 
structure of CeCoIns. So we consider that the same mechanism can be applied to the mecha- 
nism of the superconductivity of pressured CeRhIns(P > I.6GPa). 

In conclusion, we have reformulated the Eliashberg's equation for the superconducting 
transition within the quasi-particle description which takes into account of the heavy elec- 
tron mass. We also presented a microscopic mechanism for the d-wave superconductivity in 
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Celr^Coi-xIns and attempted to explain the x-dependence of T c , by using such a renormalized 
formula and the perturbation theory. 
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aG 




=a(a 3 CGGG + --) 
= (a 2 r c) 2 GGG +-■ 

Fig. 1. The renormalizcd self-energy S(k, uj) can be re-written as the expansion with respect to the 
quasi-particle Green's function G and the quasi-particle s-wave scattering interaction ri oc = a 2 ri oc 
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r(l,2;3,4)=a 2 (r(l,2;3,4)-ry 
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l u 

=a 2 ( aoQaG + 



rfoc 

=a | aG 

r]oc 

=a 2 (a 2 CGG + a 4 riGGGG + ---| 
=(aT7 0C ) GG + (aT7 0C ) GGGG + - - - 

Fig. 2. The full vertex T(l, 2; 3, 4) can be re- written as the expansion with respect to the quasi-particle 
Green's function G and the quasi-particle s-wave scattering interaction ri oc = a 2 ri oc 






Fig. 3. The Feynman diagrams of the normal self-energy up to third order. The solid and dashed 
lines correspond to the bare Green's function and the interaction, respectively. 
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Fig. 4. The Feynman diagrams of the effective interaction up to third order. The solid and dashed 
lines correspond to the bare Green's function and the interaction, respectively. The diagrams 
enclosed by the dashed line are vertex corrections. The other diagrams are included in RPA. 
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Fig. 5. The calculated T c as U is varied and for various values of n as shown in the figure. TOPT 
and TOPT(RPA-like only) in this figure mean that full diagrams and only RPA-like diagrams of 
anomalous self-energies up to third order are included, respectively. 
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Fig. 7. The density of states as n is varied, at U = 5.5 and T = 0.02. The inset shows the details 
near the Fermi level. 
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Fig. 8. The real part of the normal self-energy at the Fermi momentum, at n = 0.365, T = 0.02 and 
for various values of U as shown in the figure. 
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Fig. 9. The imaginary part of the normal self-energy at the Fermi momentum, at n = 0.365, T = 0.02 
and for various values of U as shown in the figure. 
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